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Abstract 

We consider a recursive scheme for defining the coefficients in the operator 
product expansion (OPE) of an arbitrary number of composite operators in the 
context of perturbative, Euclidean quantum field theory in four dimensions. Our 
iterative scheme is consistent with previous definitions of OPE coefficients via the 
flow equation method, or methods based on Feynman diagrams, ft allows us to 
prove that a strong version of the “associativity” condition holds for the OPE to 
arbitrary orders in perturbation theory. Such a condition was previously proposed 
in an axiomatic setting in [1] and has interesting conceptual consequences: 1) 
One can characterise perturbations of quantum field theories abstractly in a sort 
of “Hochschild-like” cohomology setting, 2) one can prove a “coherence theorem” 
analogous to that in an ordinary algebra: The OPE coefficients for a product of 
two composite operators uniquely determine those for n composite operators. 
We concretely prove our main results for the Euclidean quantum field theory, 
covering also the massless case. Our methods are rather general, however, and 
would also apply to other, more involved, theories such as Yang-Mills theories. 


1 Introduction 

There exist many different approaches to quantum field theory. Many of these attempt 
to isolate within quantum field theory a kind of algebraic skeleton, which, in a sense 
depending on the particular framework, defines the theory and dictates its properties. 
The earliest manifestation of this kind of framework is that of local quantum physics 
due to Haag and Kastler [2] which is based on nets of local algebras of operators. A 
framework to isolate the algebraic core of many 2-dimensional conformal field theories 
is the theory of vertex operator algebras [3, 4], The main idea of this framework is to 
formalise the properties of the operator product expansion (OPE) in such theories in 
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order to build an algebraic structure capable of describing many interesting models in 
two dimensions. 

Since the OPE ought to exist in any local quantum field theory in any dimension [5], 
it seems reasonable to define a quantum field theory by it, or more precisely, to attempt 
to build a self-consistent algebraic structure out of the OPE that can define a quantum 
field theory. The OPE is the statement that given a complete set of local operators 0 A , 
and given any sufficiently well-behaved quantum state X P, one has 

(OaM) ' ' ' ^A JV (*v))'F ~ ^ C\ x An {x 1 , . . . , X N ) (0 B (Xn))'Y ■ (1-1) 

B 


Here, C Ai A , are functions (or rather distributions), called OPE coefficients, and the 
symbol “ ~ ” indicates that the relation is expected to hold asymptotically at short 
distances, in the sense that the difference between the left and right hand side of (1.1) 
vanishes if %,■ —> Xn for all i < N. In models of perturbative quantum field theory, such 
as the Euclidean ^-theory, the OPE was found to be not only asymptotic, but even 
convergent, in the sense that the sum over B in (1.1) converges even for any finite 
separation of X\,..., x N [6, 7]. 

These results strongly suggest that it should indeed be possible to view the OPE 
coefficients as defining the algebraic skeleton of the theory, and the 1-point functions 
{Ob{x Ar))vp as carrying all the information about the state. The theory, then, should be 
defined by the OPE coefficients, whereas specific physical setups should be described 
by the collection of all 1-point functions, much in the way as a classical field theory 
is defined by a partial differential equation, and specific physical setups are described 
by boundary- or initial conditions for determining a given solution. (As an aside, let 
us point out that this viewpoint is, in fact, not only remarkably close to standard 
applications of the OPE in deep inelastic scattering, but also very attractive in curved 
spacetimes [8, 9], because it is much less clear there what physically preferred states 
would be in general.) 

Of course, in order to define a concrete field theory, one must have a way to 
determine the OPE coefficients in the first place. The traditional way in Lagrangian 
field theory is to go back to correlation functions and proceed e.g. by the well-known 
(perturbative) methods described in [10, 11]. This is not really satisfactory if one wants, 
as we do, to view the OPE coefficients as the primary objects defining the theory, and 
not Lagrangians or correlation functions. In order to get around this, one clearly needs 
extra information on the OPE coefficients. One central property (formalised e.g. in 
the setting [1]) is a kind of associativity (also called “factorisation" or “consistency") 
condition, which can be motivated heuristically as follows: Consider an operator 
product <9a 1 (vi)<9a 2 (v2)<9 / i 3 (^3), where x t e R 4 , and assume that x 2 is closer to X \ than 


to X 3 , i.e. 


1*1 ~ *21 

1*2 - * 31 


< 1 . 


( 1 . 2 ) 


Since the OPE is by its very nature a short distance expansion, one may hope to be 
able to perform the OPE of only the product 0 Al (x t )0 An (x 2 ) around the point x 2 first, 
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leaving 0 Ai (x 3 ) as a “spectator". Such an expansion would have the form 


<' 0 Ai (Xi)0 A2 (x 2 )0a 3 (X3 )) ~ ^C AiA2 (x 1 ,X 2 )(Ob(X2)Oa 3 (X3)) 

B 

~ '^ J C AlA2 (Xi, X 2 ) Cb A 3 (X2, X3)(Oc(x?,)) , 
B,C 


where we performed a second OPE in the second line. Comparison with eq.(l.l) yields 
an associativity condition 

C Ai a 2 a 3 (*i>*2,*3) — Ai a 2 (x\,X2)Cca 3 (X2i X3). (E4) 

c 


This condition puts strong restrictions on the OPE coefficients of the theory. To see 
this, assume also that 


1*2 ~ -pi 

1*1 - *31 


< 1. 


(1.5) 


We can repeat the argument above and arrive at the relation 


^AiA 2 A 3 (*1’ *2’ * 3) /* ] C a 2 A 3 (*2> * 3) ^ Al c(^l’ *3) • 

c 


( 1 . 6 ) 


The requirement of consistency of the alternative expansion schemes (1.4) and (1.6) on 
the domain 0 < \xi - X 2 I < 1*2 - * 3 ! < |*i - * 3 ! yields 

Yj C MaM i» * 2 ) C®a 3 (* 2 , * 3 ) = Yj C a 2 a 3 (*2, * 3 ) C AlC (x 1 , *3) , (1.7) 

c c 

which encodes highly non-trivial relations between the OPE coefficients. It was shown 
in [1] that these have various consequences: 

• Multipoint OPE coefficients C A A v are uniquely determined in terms of the 
two-point coefficients C A]A2 - 

• Deformations (=perturbations) of OPE coefficients can be characterised as a 
cohomology of Hochschild type. 

• OPE coefficients can be viewed as a (non-conformal, higher dimensional) version 
of vertex operator algebras. 

The formal “derivation" of the associativity condition presented above is, of course, far 
from rigorous: For one thing, we have introduced the OPE as an asymptotic expansion, 
but in (1.2) and (1.5) we demanded finite separation of the points X\,X 2 , *3. Furthermore, 
it is not obvious in what sense, if at all, the partial OPE performed in (1.3) holds. Lastly, 
we have implicitly exchanged the order of two infinite series in the step from (1.3) to 
(1.4) without any justification. Nevertheless, it is possible to see in some non-trivial 
examples of field theories such as in the massless Thirring model [12], or in the context 
of 2 dimensional conformal field theories [13] that the strong form of the associativity 
condition (1.7) in fact holds. Unfortunately, the arguments presented in these works 
are very specific to the peculiar properties of such models, giving no hint whatsoever 
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what the situation might be e.g. for perturbatively defined models in Lagrangian field 
theory. 

In the present paper we show that associativity of the OPE indeed holds to all orders 
in the perturbative Euclidean ^-theory. In fact, we even prove a generalisation of 
eq.(1.4) to more than three fields: 

Theorem 1. Denote by [A] the dimension of the composite field Oa- At any perturbation 
order r 6 N in Euclidean theory ; there exist constants c, K > 0 such that 




.,x N )~ ^ C c Ai .._ Am {xu---,Xm) CcA m+1 ..A n ( X M, Xm+\, . . . , Xn) 
[C\<D 

<Kf 


r-th order 


D ±1 
2 


n , 0 nc-(2 ,[A, ]+m max(^, \x t - va,|) [B]+1 

D + Z \ \<i<N m 


min | X; - w|2/ A d +1 

l<i<j<N 


( 1 . 8 ) 


holds for any x\,..., x^ such that 


maxi<,< M \xj - x M \ 

m i n M<j<N \Xj — x M I 


(1.9) 


where c = c(r ) and K = K(r, A i,... ,A n , B ) do not depend on D. Since the r.h.s. of (1.8) 
vanishes in the limit D —> oo, the bound implies that the associativity property 

Cai..a n ( x 1’ ■ ■ ■ ’ x n) = ^ • ■ ■ > x m)C CA m+1 An (xm, ■ • • ,Xn) ( 1 - 10 ) 

[C]<D 


holds up to any perturbation order on the domain defined by (1.9). 


Remark: A much weaker version of associativity was previously derived in [14], 
There, it was shown that eq.(l.lO) indeed holds up to any perturbation order, but only 
on the smaller domain 

max,< ; < M \x t - x M \ 

0 c— : --- r <£ (Eli) 

min />M I Xj - x M \ 

for some constant 0 < e c 1 which moreover decreases with the perturbation order. 
The weaker version is not suited in order to derive (1.7). Furthermore, the weaker 
version gives the misleading impression that associativity breaks down altogether 
beyond perturbation theory. 

This result suggests that a quantum field theory can be defined by a set of OPE coeffi¬ 
cients satisfying (1.10) on the domain (1.9), together with other simple straightforward, 
and reasonable requirements, see section 2 (for more details see [1] and also [15, 16] 
for curved spacetimes). 

Even though, thanks to the above theorem, we may now feel much more confident 
that this viewpoint on QFT is correct, it does not tell us how to actually find QFTs, 
i.e., how to find actual solutions to the consistency requirements (1.10). Here a further 
independent idea is needed. This idea is to investigate how, given one solution to the 
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consistency relations (e.g. the Gaussian free field), one can deform this solution to 
another one. As we recall below, one can nicely formulate an abstract deformation 
(=perturbation) theory of the algebraic structure based on (1.10) wherein perturbations 
are characterised as elements of some Hochschild type cohomology ring. However, 
this still does not give a good practical way of actually finding perturbations (to all 
orders in some small parameter, or even finite ones). Instead, we are going to rely on a 
recently found recursion formula for perturbative OPE coefficients [17]. This recursion 
formula is derived from the differential equation (a caret" denotes omission) 

dgC B Al ...A N (x u ...,x N ) = -J d 4 y [cf Al A)V (y, xi,, x N ) 

N _ 

- Yj C ZAi(y’ X ^ C A l ...A i D...A N ( ' XU C M-A N ( X l’ 

i= 1 tD]<[A;] [D]<[B] 

for the change of an OPE coefficient if we change the action of the theory by a term of 
the form gOg (where gO A would be gtp 4 in our model). It is this relation, together with 
the well-known formulae for the OPE coefficients of the free theory (g = 0), which 
is used in this paper to construct the coefficients of the interacting theory order by 
order in g, and to prove theorem 1. The bottom line is that this recursion formula (or the 
differential equation), together with the consistency relation (1.10) completely determine 
the OPE coefficients of a theory - hence the theory itself - and that these conditions are 
mutually consistent with each other. 

This paper is organised as follows: We put our results into the context of axiomatic 
approaches in section 2. Section 3 contains the main results of the paper, which are 
then proved for the case of massive fields in section 4. The generalisation of the proof 
to massless fields can be found in section 5, followed by our conclusions in section 6. 
Some technical estimates are moved to an appendix. 


• • • i XNiCgpiy, yv)J > 

( 1 . 12 ) 


2 General framework for QFT and remarks 

Before delving into the derivation of the main results of this paper, we would like to 
explain the wider context provided by a specific proposal for the structure of QFT [1]. 


OPE algebras: This framework is intended to formalise the properties of the OPE. 
In order to avoid writing many indices, one associates local fields 0 A in the theory 
with vectors |va) in some abstract vector space called V. The space V is assumed to be 
graded in various ways which reflect the possibility to classify the different composite 
quantum fields in the theory by their spin, dimension, Bose/Fermi character, dimension 
etc. Thus, for example, if Vd is the space of all fields of a fixed dimension D, then 




D 


Vd- 


(2.13) 


The infinite sum in this decomposition is understood without any closure taken. In 
other words, a vector |v) in V has only non-zero components in a finite number of the 
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direct summands in the decomposition (2.13). Typically the set of possible Z)-values is 
discrete and each dim Vd < 00 1 - 

On the vector space V, we assume the existence of an anti-linear, involutive opera¬ 
tion called ★ : V —> V which should be thought of as taking the hermitian adjoint of the 
quantum fields. We also assume the existence of a linear grading map y : V —> V with 
the property y 2 = id. The vectors corresponding to eigenvalue +1 are to be thought 
of as "bosonic", while those corresponding to eigenvalue -1 are to be thought of as 
"fermionic". 

So far, we have only defined a list of objects—in fact a linear space—that we think of 
as labelling the various composite quantum fields of the theory. The dynamical content 
and quantum nature of the given theory is next incorporated in the OPE associated 
with the quantum fields. This is a hierarchy denoted 

C = (c(- -), C(- -), C(- ( 2 . 14 ) 

where each (x \,..., x N ) h-> C(x 1 ,, x N ) is a function on the "configuration space" 

M n := {(vj,..., x N ) € (R 4 ) w I Xi Xj for all 1 < i < j < N) , (2.15) 

taking values in the linear maps 


C(xi,... ,x N ) : V® • • • 0 V —> V, (2.16) 

where there are N tensor factors of V. (The range of C(x \,..., x N ) is actually in the 
closure V** of V but we do not distinguish this in our notation.) The components of 
these maps in a basis of V correspond to the OPE coefficients mentioned in the previous 
section. For one point, we set C(x 1 ) = id : V —» V, where id is the identity map. 

In order to have any chance of imposing stringent consistency conditions of the 
nature described in section 1, the maps C(—,...,-) must be real analytic functions on 
M n , in the sense that their components Cf An ( x i> • • • > x n) := (vb\C(x \,..., Xn)\va { ® 

... ® Va n ) are ordinary real analytic functions on with values in C. The basic prop¬ 
erties of quantum field theory are then expressed as the following further conditions 
on the OPE coefficients: 


C1) Hermitian conjugation: Denoting by ★ : V — » V the anti-linear map given by 
the star operation, we have [4r,y] = 0 and 

C(x = *C(x i,...,x N ) (2.17) 

where -k® N := ★ ® • • • ® ★ is the Wfold tensor product of the map ★, and where 7 
denotes complex conjugation. 


fin order to have a reasonable theory possessing sufficiently many states it is natural to demand a 
finiteness property of the kind ° dim Vq < oa for 0 < q < 1. 
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C2) Euclidean invariance: For a suitable representation R of Spin(4) on V and 
a G EL 4 , g G Spin(4), we require 

C(gx i + a,.. ,,gx N + a) = R*(g)C(x u ..., x N ) R(gf N , (2.18) 

where R(g)® N stands for the A-fold tensor product F(g) ® • • • ® F(g). 

C3) Bosonic nature: The OPE-coefficients are themselves "bosonic" in the sense 
that 

C(x u ...,x N ) = yC(x l ,...,x N )y® N (2.19) 

where y® N is again a shorthand for the w-fold tensor product y ® • • • ® y. 


C4) (Anti-)symmetry: Let r,_i ; = (i — 1 i) be the permutation exchanging the 
(i - l)-th and the j-th object, which we define to act on V ® • • • ® V by exchanging the 
corresponding tensor factors. Then we have 

C(x 1,..., JC/-1, x u ..., Xjv) T i- 1 ,/ = C(x 1,..., X/, ..., x N ) (2.20) 

F, := ^ id m ~ X) ® (id - y) ® id* (A/ - /} . (2.21) 

for all 1 < i < N. Here, the last factor is designed so that bosonic fields have symmetric 
OPE coefficients, and fermionic fields have anti-symmetric OPE-coefficients. The last 
point Xn and the A-th tensor factor in V ® • • • ® V do not behave in the same way 
under permutations, and the formula has to be slightly altered. See [1, eq.(3.38)] for the 
corresponding formula. 


C5) Scaling: Let dim : V —> V be the “dimension counting operator”, defined to act 
by multiplication with D G BL+ in each of the subspaces V D in the decomposition (2.13) 
of V or, put differently, dim v A ) = [A] • |va). Then we require that 1 G V is the unique 
element up to rescaling with dimension dim(l) = 0, and that [dim, y] = 0. 

Furthermore, we require that, for any 6 > 0 and any (x \,..., x„) G M n , 

lime [Al]+ - +[Aw] - [B]+,5 C^... AiV (exi,..., ex N ) = 0. (2.22) 


C6) Identity element: We postulate that there exists a unique element 1 of V of 
dimension [1] = 0, with the properties 1* = 1 , y( 1) = 1 , such that 


C{ x u ...,x N )\v l ® • • • 1 ® • •• Vjv_i> = C(x x h ...x N )\v { ® • • • ®Vjv_i) . (2.23) 

where 1 is in the i-th tensor position, with i < N — 1. When 1 is in the iV-th tensor 
position, the analogous requirement takes a slightly more complicated form (see [1, 
chapter 3]). 
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C7) Factorisation: 

C(xx N ) = C(x M ,..., x N )(C(x\ ,...,x M )® id® XN ~ M) ) (2.24) 

on the domain 

maxi<,< M \xj - x M \ < l ^ 

mi n M<j<N \ x j ~ x m\ 

Note that this condition is an “index free" restatement of (1.10), the main result of our 
paper in the context of perturbation theory. 

Definition 1 . A quantum field theory is defined as a pair consisting of an infinite dimen¬ 
sional vector space V with decomposition (2.13) and maps 4r,y, dim with the properties de¬ 
scribed above, together with a hierarchy of OPE coefficients C := (C(-, -),C(-, -),...) 

satisfying properties Cl)-C7). 

It is natural to identify quantum field theories if they only differ by a redefinition 
of the fields. Informally, a field redefinition means that one changes ones definition 
of the quantum fields of the theory from O a (x) to 0 A (x) = £ B Z^<9 e (v), where Zf is 
some matrix on field space. The OPE coefficients of the redefined fields differ from the 
original ones accordingly by factors of this matrix. We formalise this in the following 
definition: 

Definition 2. Let (V, C) and (V, C) be two quantum field theories. If there exists an 
invertible linear map Z : V —> V with the properties 

Z R(g) = R(g) Z, Zy = yZ, Z+ = *Z, Z( 1) = T, dimZ>Zdim, (2.26) 

together with 

C(x x ,...,x n ) = Z- x C(x u ...,x n )Z® n (2.27) 

for all N, where Z m = Z ® • • • ® Z 3 then the two quantum field theories are said to be 
equivalent, andZ is said to be afield redefinition. 

The main result of our paper, Thm. 1, leads to the following conclusion: 

Corollary 1 . The OPE in perturbative Euclidean ipff theory satisfies axioms Cl)-C7) in 
the sense of formal perturbation series in g, i.e. at each fixed order in g. 

Proof. The symmetry requirements Cl)-C4) and the identity axiom C6) are quite easily 
checked: They can be explicitly checked in the free theory, and one verifies directly that 
they are preserved by the recursion formula (1.12), which we use to define perturbative 
OPE coefficients. The scaling requirement C5) follows e.g. from the bounds proven 
in [17]. By far the most non-trivial challenge is to prove C7) (factorisation). This is the 
content of thm.l of the present paper. □ 
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Vertex algebras: Another corollary of theorem 1 is that perturbation theory defines 
an analog of a vertex operator algebra: First, define vertex operators Y(x, v) : V —> V as 
the endomorphism of V whose matrix elements are given by 

(v c \Y{x, v A )\v B ) := C AB (x , 0) (2.28) 

for any x ± 0. The relation (1.7), which is a consequence of our main theorem, may 
now be written as 

Y(v,x)Y(w,y ) = Y(Y(v,x-y)w,y), (2.29) 

where the spacetime arguments are required to satisfy |jt| > |y| > |x - y| > 0 and where 
v, w are elements of V. An almost identical quadratic relation first appeared in the study 
of conformal field theories in two dimensions, where it is one of the crucial properties 
(called “locality condition") of the vertex operator algebras [4], It should be stressed, 
however, that in our context, where conformal symmetry is not required, the condition 
above is really a highly non-trivial statement on the convergence of the infinite sums 
implicit in eq.(2.29), whereas the same equality in the CFT context is understood in 
terms of formal power series. 

Abstract perturbation theory: The constraint imposed by the factorisation condi¬ 
tion C5) at the three point level can be rewritten as 

C(x 2 , x 3 )[C(x u x 2 ) <8> id ) - C(x u x 3 )[id <g> C(x 2 , x 3 )) = 0 

for 0 < |a'i - x 2 \ < \x 2 - * 3 1 < |xi - * 3 1, 


which is just an “index free" version of eq.(1.7). Although we will not use this in the 
present paper, all higher constraints can be derived from this one, see [1], In the 
very abstract general framework of an OPE algebra, we may ask the question when 
it is possible to find a 1-parameter deformation C{x\, x 2 ;g) of these coefficients by a 
parameter g so that the associativity condition continues to hold, at least in the sense 
of formal power series in g. (Actually, the analogues of the symmetry condition (2.20), 
the scaling condition (2.22), the hermitian conjugation, the Euclidean invariance, and 
the unit axiom should hold as well for the perturbation. However, these conditions are 
much more trivial in nature than (2.30), because the conditions are linear in C(xi,x 2 ). 
These conditions could therefore easily be included in our discussion, but would distract 
from the main point.) 

One can show that such perturbations can be characterised in a cohomological 
framework. To set up this framework, we consider the non-empty, open domains of 
(R 4 )^ defined by 

Fn = {(-*T • • • > -Tv) 6 M n ', Hi-i < r i-u < ri_ 2i < ■ • • < r 3i , l<i<N}cM N , (2.31) 

where r,y := \x t - Xj\. We define Q N (V) to be the set of all real analytic functions / v on 
the domain T'n that are valued in the linear maps 

In(x i ,..., x N ) : V ® • • • 0 V -> V, (x u ..., x N ) 6 T N ■ (2.32) 
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We next introduce a boundary operator b : Q N (V) —> Q /v+1 (V) by the formula 


(bf N )(x u ..., x N+l ) := C 0 (*i, x N+ i)(id 0 / w (x 2 ,..., *az+i)) 

N 

+ ^(-1)7^1, ® C()(x,-, X,+ i) 0 i(f (N ~ l) ) 

i=l 

+( _ 1) /V+1 CoCw* Tv+i)(./v7i > • • •, -Tv) ® id) . (2.33) 


Here Co(-H> * 2 ) is the OPE-coefficient of the undeformed theory defined by g = 0, and 
a caret means omission. The definition of b involves a composition of Co with fy, and 
hence, when expressed in a basis of V, implicitly involves an infinite summation over 
the basis elements of V. We must therefore assume here (and in similar formulas in the 
following) that these sums converge on the set of points (x\,..., Xn+i) in the domain 
T'n+x- We shall then say that bf N exists, and we collect such f N in the domain of b, 

dom(b) = € f2 w (y) | bfN exists and is in f2 A,+1 (V)}. (2.34) 

N>1 

When we write bf N , it is understood that f N 6 Q v ( V) is in the domain of b. One can 
show: 


Lemma 1. The map b is a differential, i.e., b 2 f N = 0 for f N in the domain ofb such that 
bfN is also in the domain ofb. 

Let us define the kernel Z N (V,C) of b on Q N (V) as the linear space of all 6 
Q N (V) n dom(b) such that bfN = 0. Similarly, define the range B N (V,C) in Q N (V) to be 
the linear space of all f N = b/'v- 1 such that f N ~ i 6 Q V_I (V) n dom(b) and such that f N 
is in dom(b). By the above lemma, we can then define a cohomology ring associated 
with the differential b as 


Z N {V\C) _ {kerb : Q. N (V) -> D w+1 (k)} n dom(b) 
B N (V; C ) {ranb : 0^(7) -> Q W (V)} n dom(b) ' 


As we will now see, the problem of finding a 1-parameter family of perturbations 
C(x i, jc 2 ; g ) such that our associativity condition (2.30) continues to hold for C(x i, jc 2 ; g ) 
to all orders in g can be elegantly and compactly formulated in terms of this ring. If we 
let 

1 d‘ 

Cfx,,x 2 ) = - —C(x l ,x 2 ;g) , (2.36) 

z! dg‘ 


8=0 


then we note that the first order associativity condition, 


C 0 (x 2 , x^Cxixuxf) 0 id) - C 0 (x u x 3 )[id 0 Ci(.r 2 , x 3 )J + 
Ci(x2,x 3 ){Co(xux 2 )®id) -Ci(xuX3)(id®Co(x2,x 3 )) = 0 , (2.37) 


valid for {x\,x 2 , xf) 6 T 3 , is equivalent to the statement that 


bC i=0, 


(2.38) 
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where here and in the following, b is defined in terms of the unperturbed OPE- 
coefficient Co- Thus, C\ has to be an element of Z 2 (V\Cq). Let Z(g) : V —» V be 
a g-dependent field redefinition in the sense of defn. 2, and suppose that C(x\, x 2 ) and 
C(x i,x 2 ; g) are connected by the field redefinition. To first order, this means that 

Ci(xi,x 2 ) = ~Z\Cq(X] , x 2 ) + Co(xi, x 2 )(Zi 0 id + id 0 Z 1 ), (2.39) 

or equivalently, that hZ\ = C\, where Z, = j-;Z(g)\ g= o. Thus, the first order defor¬ 

mations of Co modulo the trivial ones defined by eq. (2.39) are given by the classes in 
H 2 (V',Co). The associativity condition for the z'-th order perturbation (assuming that 
all perturbations up to order i - 1 exist) can be written as the following condition for 

(x u x 2 ,x 3 ) e Ty. 

Co(* 2 , x 3 )(Cj(xi, x 2 ) 0> id} - Cj(x 1 , x 3 ){id 0 C 0 (x 2 , x 3 )} + (2.40) 

Cj(x 2 , v 3 )(c 0 (vi, v 2 ) 0 id} - Co(xi,x 3 )(id 0 Cfix 2 , x 3 )J = wfx u x 2 , x 3 ), 

where w t 6 fl 3 (V) is defined by 

(-1 

Wi(x 1 , x 2 , x 3 ) := Ci-j(x 1 , x 3 )(id 0 Cj(x 2 , x 3 )) - Ci-j(x 2 , x 3 )(C j (x\ , x 2 ) 0 id). (2.41) 

i=i 

We assume here that all infinite sums implicit in this expression converge on ( F 3 . This 
equation may be written alternatively as 

bCi = wi. (2.42) 

We would like to define the z-th order perturbation by solving this linear equation 
for C;. Clearly, a necessary condition for there to exist a solution is that bw, = 0 
or Wi 6 Z 3 (V, Co), and this can indeed be shown to be the case. If a solution to 
eq. (2.42) exists, i.e. if Wj e B 3 (V,Cq), then any other solution will differ from this 
one by a solution to the corresponding "homogeneous" equation. Trivial solutions 
to the homogeneous equation of the form bZ, again correspond to an z'-th order field 
redefinition and are not to be counted as genuine perturbations. In summary, the 
perturbation series can be continued at z'-th order if [w ( ] is the trivial class in H 3 (V;Cq), 
so [Wj] represents a potential z'-th order obstruction to continue the perturbation series. 
If there is no obstruction, then the space of non-trivial z'-th order perturbations is given 
by H 2 (V',Co). In particular, if we knew e.g. that H 2 {V\Cq) 4- 0 while H 3 (V',Co) = 0, 
then perturbations could be defined to arbitrary orders in g. 

The relationship of this abstract framework with the results of the present paper is 
the following: 

Corollary 2. Let Co be the OPE coefficients of a free, scalar Euclidean quantum field 
theory, and letCj, j > 0 be their perturbations, as defined by the recursion formula (1.12). 
Then 

a) Ci is a non-trivial element o/7/ 2 (V;Co), and 

b) all higher obstructions [wi] £ H 3 (V;Co) vanish. 
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Proof. Non-triviality of Ci follows from the fact that the recursion formula (1.12) can 
not be written as a mere redefinition of the composite fields. The second point, i.e. 
vanishing of obstructions [w ( ] 6 H 3 (V;Co), follows directly from the main result of the 
present paper, thm.l, because it is equivalent to associativity order-by-order in g. □ 

3 The Associativity Theorem 

In the present section we are going to state our other main results, which will imply 
the bound stated in thm. 1 within perturbative Euclidean ^-theory in four dimensions 
with classical action 

S = Jd 4 x^(d(f) 2 + + |^ 4 J . (3.43) 

Throughout the present section we will restrict attention to the massive case m 2 > 0. 

The generalisation of our proof to massless fields is discussed afterwards in section 5. 
We write the composite operators of our model explicitly as 

0 A =d a 'tp---d a "ip, A = (a u ...,a„), ctj e N 4 , (3.44) 

which means that the corresponding dimension of the field Oa is given by 

n 4 

[A] = ^(1 + |ari|), where \a\ = ^ \af\ for aeN 4 . (3.45) 

i=\ fi=l 

Let us denote the (formal) perturbation series for OPE coefficients by 

oo 

C B Al ...A N ^x, ...,x n )=: YjCr)M...A N (.x 1’ • • • > x N ) ■ g r , (3.46) 

r =0 

where the perturbative OPE coefficients (C r ) Al An are defined recursively through 
eq.(1.12). Further, denote by 

(•‘ ’ T/v) • — 

(Cr)A l ...A N ( X l’---’ X N)~ XX (C,4. ..AmO*-!’ • • ■ ’ x m) (Pt) C A u+i ... A n ( X M, ■ ■ ■ > Xn) 

s+t=r [C\<D 

(3.47) 

the remainder of the associativity condition at r-th perturbation order and truncated 
at operators (9 C of dimension [C] = D. Our strategy is to establish the bound (1.8) by 
an induction which is based on the recursion formula (1.12). In order to obtain the 
sharp bound (1.8), we will have to formulate our induction hypothesis not in terms of 
the remainder functions (R?) Al a m -a m +i a n > terms of much more general objects, 

containing multiple summations over products of OPE coefficients (see definition 4 
below). These more general expressions are most conveniently organised in terms of 
decorated rooted trees. Before we can state our main inductive bound, we therefore 
have to introduce some additional notation. 

First, we agree on a vocabulary for rooted trees T, which is summarised in the 
following glossary (cf. [18, chapter 3.2.2]): 
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Symbol Definition 


'V(T) Vertices of the tree T. 

TXT) Leaves of T, i.e. vertices of degree 1 (the degree of a vertex is the number 

of edges adjacent to it). 

K(T) The root of T,<R £ 'V. 

I(T) Internal vertices of T, i.e. non-leaf vertices. 

Ik(T) Internal vertices of T without the root, i.e. I<r := I \ % . 

B{T) The set of branches of T. A branch b £ B is a path connecting a leaf to 

the root, where we use the convention that leaves and root are not part 
of the branch, i.e. Be. 1%. 

ch(v) The children of a vertex veT are the vertices adjacent to v which are 

further from the root. 

pa(v) The parent of a vertex v 6 < V is the vertex adjacent to v which is closer 

to the root. 

sb(v) The siblings of a vertex v £ f V are the children of the parent of v not 

including v itself, i.e. sb(v) := ch(pa(v)) \ {v}. 
de(v) The descendents of a vertex veT are the vertices on the paths from v 

to the leaves. 

an(v) The ancestors of a vertex veT are the vertices on the path from v to 

the root. 

Next, we add decorations to these trees: 

Definition 3 (Weighted trees). Let x = (x\, ... ,x^) £ M. 4N and A = (Aj,... ,A„), where 
Aj 6 N 4 "' are multi-indices and where n > N. We define T{x\ A) to be the set of rooted 
trees T with the following properties: 

1. T has n vertices and N leaves. 

2. Vertices in I<n have degree larger than 2. 

3. To each vertex v £ < V(T) we associate a pair (x v , A v ) called the weight ofv, where 
x v £ R 4 is a four-vector andA v £ N 4 "' a multi index, such that 

• ifv £ I(T), then x v £ \x w : w £ ch(v)}, i.e. x v has to be equal to one of the 
four-vectors associated to the children ofv. To the leaves v e -C(T) we associate 
bijectively the vectors (x\,..., Xn), i.e. (x v ) ve £ = x. 

• {A v ) ve y( T ) = A, i.e. the mapping between multi-indices and vertices is one-to- 
one. 

See fig.l for an example of three such trees. 

With this notation in place, we can now give a compact definition of the objects 
appearing in our induction hypothesis: 

Definition 4 (Contractions of OPE coefficients). Given a tree T gT(x; A), we define 
(P, )(T) := J - 1 (C r f ( f e)eeMv) ((x e ) e€C h( V )', T>) • (3.48) 

vel(T) 2 r u =r 
ueI(T) 
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(B,x n) 





Figure 1: Example of weighted trees To e T((A\, ..., Ajq, B)\ (x \,..., x^)), T\ 6 
T{{A U ---,A N .B,C)\{x\,...,x N )) and T 2 eT((Ai,...,A 5 ,B,Ci,C 2 );(xi,...,x 5 )). 


The argument behind the semicolon in the OPE coefficients specifies the reference point, i.e. 

Ca 1 ...a jv (- x ' 1 ’ • ■ • ’ x n, x \) = Qa 2 ... AnAi ( x 2, • ■ ■, Xn, X\) (3.49) 


for example. 


Examples: For the weighted trees depicted in fig.l, the definition yields 

(Prm) = <Pr'fi lA2mAK (x 1 ,X 2 ,... ,x N ) (3.50) 

(Pr)(T\ ) = ^ X m) (Cr 2 ) C A M+1 ...A N ( X M, ■ ■ ■ , Xn) (3.51) 

r\ +r 2 =r 

(Pr)(T 2 ) = iC r p ( - A \ A fx\,x 2 ) (CrjfiffiM, X 5 ) (C ) . 3 )c 1 A 3 C 2 fe, *3, x 5 ) . (3.52) 

ri+r 2 +r 3 =r 


We are now ready to state our second main theorem, which directly implies theorem 1. 


Theorem 2. Up to any perturbation order r 6 N, OPE coefficients of massive Euclidean 
(p^-theory satisfy the following two properties: 

(a) Given a tree T e T{x\ A) and given a collection of integers (D,) /e j K , fix any branch 
b 6 S(T) in the tree such thaf x v = x w for all v, w 6 b(T) and define the shorthand 
D t := Ya V ££\jr(T)\Av\- For any choice of constants s e (0,2~ (I)r+4r+3) ] andS v 6 (0,1), 
one has the bound 


n z 

ieI K (T) [A,]=A- 


max | Xi - x%\ [A k] 

iech(ft) 


[1 min |Xi - Xj\ [A - ] 

ie£ t' esb 9') 


I h " 1 * 

ieln 


x K 


(1 + sffi" 


b[A w ] 


l''Xvc'V\b\Av ] 


• Y](D W + lf r ) 8 

web 


><n 


SU P/«h( V )(l*i - *vl, 1M) 

m e(A > ) min |.y - r.| 1+ ^ A >) 

V i*fe ch(v) 


\<5v 


(3.53) 


2 Such a branch exists for every tree T. In fact, it is not hard to see that the number of such branches 

is equal to the degree of the root vertex 'fi of T. 
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where the constant K > 0 depends neither on the integers D t nor on s or the 6 V , 
where 


max 


eech(/) 


l-T ~ Xj | 


UT) 

mm e6sb(i) \x e - Xi\ 

where 6 is the Heaviside step function 3 and where A,, := 2 eech (v)[A e ] - [A,,]. 
(b) For any £ < 1 one has 


(3.54) 


...,X N ) = 0, 


Z)—> oo 


(3.55) 


where £ is defined as in eq.( 1.11). 


Remark: Before we come to the proof of the theorem, let us take a moment to 
have a closer look at the result in order to get a better intuition for the complicated 
expression (3.53). The origin of the various terms in the bound (3.53) can be roughly 
understood as follows: 

1. The first line reflects the behaviour one would expect from naive power counting 
if one assumes that an OPE coefficient behaves as Aw (t. • ■ •, -Tv) ~ max |t - 

TvP 1 / FL m i n It - t/F [A,] . 

2. The second line captures all combinatorial factors, in particular those caused by 
the summations over multi indices [C,] = D, associated to the internal vertices 
of the tree T and those arising in perturbation theory. Note that only this second 
line depends on the perturbation order r. 

3. In the last line, the factors including the Heaviside function are a relict of the 
exponential decay of the massive propagator. These factors are needed in the 
induction in order to avoid infrared divergences. Finally, the factor (sup lgch(v) (|T - 
jc v |, 1/m)/ min |t - t|) < 5, in the last line reflects the fact that naive power 

i*j e ch(v)) 

counting only holds up to logarithms once we proceed to higher orders in 
perturbation theory. We note also that the bound diverges if we set the mass m 
to zero. 


Proof of theorem 1. As mentioned in the introduction, theorem 1 can be derived straight¬ 
forwardly from theorem 2. To see this, note that eq.(3.55) implies 

)a 1 ...A m ;A m+ i...AaX- X: 1’ • • • > Tv) 

n+r 2 =r[C]>D \C\>D 

where T , 6 T({A\,... ,A N , B, C); (jci, ..., Tv)) is the tree depicted in fig.l. We can now 
use the bound (3.53) to estimate the right hand side. The infinite sum can be bounded 

3 We use the convention 0(0) = 0. 
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using the inequality 


d>D 


<(f (i +®rr 


D + 2 


1-^(1+e) 


gr+l 


8 r+1 t>j-+l 


(8' +I Dr)!, 


J3.57) 

where Dy = X|=i [A,] + [f?] and where we chose s small enough such that (1 +e) 8 " + £ < 1. 
In particular, we are free to choose for example (1 +e) 8 ’ =1/ After simple algebraic 

manipulation and absorbing some factors into the constant K, we arrive at (1.8). □ 


The reader may wonder at this stage why we derive the rather complicated bounds 
(3.53) on the objects ( P r )(T ) if we are eventually only interested in the simpler bound 
(1.8). The reason for this apparent detour lies in the fact that the bound (1.8) itself 
is not suited for the induction we are using. Roughly speaking, the main technical 
problem with an induction based on the remainder (. a m -a m+1 a n comes from the 
fact that one wants to avoid making relatively rough estimates for the summations 
over multi-indices appearing in the recursion formula (1.12). As an example, one would 
have to use estimates like 

| ^ s)A 1 ...A M ^t \*C-Am+\-A n 

[C]<Z> [C]<D 

As it turns out, such estimates lead to unwanted combinatoric factors of the form 
c D for some constant c > 1, which accumulate for every iteration of the recursion 
formula. As a result, one is led to an associativity condition that gets weaker as 
the perturbation order increases (similar to the result derived in [14], see also the 
remark below theorem 1). Our solution to this problem is to estimate the objects 
rLeJ^n Yj[Ai\=Di^Pr){T), which include multiple sums over multi-indices A,- and which 
thereby allow us to avoid weak estimates of the type (3.58), i.e. we never have to 
“pull the modulus inside the sum". The formulation in terms of rooted trees is further 
convenient in order to keep track of the various terms generated by the recursion 
formula (1.12), and in particular in order to verify cancellations of divergent terms in 
the recursion as discussed in more detail in the next section. 


< 


£ F> 


.Am 


D\B 


(K) 




(3.58) 


4 Proof of theorem 2 

In the present section we are going to present the proof of theorem 2, which proceeds 
by induction in the perturbation order r. Before we get to the details of this rather long 
line of arguments, let us give a brief overview of the general strategy and the main 
steps followed in this section. 

Induction start (sec. 4.1): Theorem 2 makes two claims, namely the bound (3.53) and 
the convergence property (3.55). Thus, our aim is to prove both these properties 
for r — 0, i.e. within the free theory. In this simple case, we can treat the problem 
explicitly using mainly Wick’s Theorem. Namely, we can write down an explicit 
representation for the zeroth order OPE coefficients [see eq.(4.60)], and we then 
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generalise this representation to the objects of interest (P,.)(T) [see lemma 2]. 
With this representation at hand, we can a) derive the claimed bounds (3.53) 
[see subsection 4.1.1] and we can b) check for convergence of the associativity 
condition [see subsection 4.1.2]. 

Induction step (sec. 4.2): Our aim is again to prove the bound (3.53) and the conver¬ 
gence property (3.55), but now at perturbation order r + 1, under the assumption 
that both these properties hold up to order r. Our main ingredient here is the 
recursion formula (1.12), which implies a corresponding recursion formula for 
the objects of interest (P r )(T) [see eq.(4.91)]. This formula allows us to establish 
bounds on (P r +i)(T) in terms of an integral over objects at order r, for which we 
can use the inductive bound by assumption [see subsection 4.2.1], In order to 
verify the bound (3.53) at order r + 1, it then remains to estimate this integral. 

Here some care has to be taken, since the individual terms under the integral 
generated by the recursion formula are in fact divergent. One has to make use 
of cancellations between such terms in the potentially dangerous integration 
regions, which can be nicely organised with the help of our tree notation. Thus, 
we decompose R 4 into various intermediate-, short- and large-distance regions, 
and we consider the integral over these regions separately. The cancellations 
between divergent terms then can be seen to follow from the associativity con¬ 
dition (3.55) at order r, and the bound (3.53) can be verified in each region by 
rather straightforward computations. 

Finally, to prove the convergence property (3.55) at order r + 1, we once again 
use the recursion formula in order to express the associativity remainder at order 
r + 1 in terms of an integral over quantities at order r. Then, using the bound 
(3.53) at order r + 1 that we have just verified, we can exchange the order of the 
integral with the limit D —» oo, which leads to a vanishing integrand, and thus to 
a vanishing remainder as claimed [see subsection 4.2.2]. 


4.1 Induction start: The free theory 

Our aim in this section is to verify the two hypotheses of theorem 2, i.e. the bound (3.53) 
and the convergence property (3.55), for free quantum fields. This will be achieved by 
giving an explicit representation for the objects (Pq)(T), which is obtained with the 
help of Wick’s Theorem. 

To derive this representation, let us start with the simplest possible trees, i.e. let 
To 6 7~(v; A) be any tree whose only internal vertex is the root (such as To in fig.l). 
Recall from our example in eq.(3.50) that the corresponding expression (Po)(To) is 
simply a single OPE coefficient. For concreteness, we write the multi indices A v e N 4 " v 
associated to the vertices v e < V(T) explicitly as, 


A v (cy 


1 >' 


,a, 


,,„ v ) ,0 Av =d a '-'<p---d a ^<p 


,a v 


(4.59) 


Wick’s Theorem then implies the convenient representation (this follows from the 
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standard definition of OPE coefficients for a free scalar field, see e.g. [17, eq. (2.56)]) 4 

(*Wo) = (C 0 )*«((x v ) v€£{To y,x mo J = V FT Z (vMwJ ) (4.60) 

^WoMKv.oJicr 

where the &v n v X n v -matrix X is defined as 


X 


(v,i)(w,j) 


y‘dZ'Mx v - Xw ) 

0 


for v,w *R,v j: w 
for v K,w = K 
for v = w, 


(4.61) 


where 


AW 


1 r e ipx 4 

(27t) 2 J p 2 + m 2 P 


(4.62) 


is the Euclidean propagator, and where s Hf('V(7o)) is the set of perfect matchings on the 
vertices {(v, O^cr^}- A perfect matching on a vertex set I is a set of edges such that 
each vertex in / is incident to exactly one edge. Figure 2 illustrates in a simple example 
how to obtain the r.h.s. of eq.(4.60) from a given tree Tq. 


rp (AftiX) 

A 


a e Tt(V(T 0 )) 

/ \ 


&V,1 &w, 1 a w ,2 Otw, 3 &7Z,1 &7Z,2 

/ \ 

{A i>,#) 


• 9 • • • • 


fua 


(otTz^i — a-w^V- 


( X -y) a K,2- a ™,3 

(&ll,2—Oi Wy 3)\ 


dx vA dy™' 1 A(x 


y ) 


Figure 2: From trees to OPE coefficients: Given a tree To [top left] with multi-index labels A v = 
(a v ,i),A w = (a Wt i,a Wt 2 , cr w ^) and A<r = (a«,i, cc<r,2), we obtain perfect matchings cr 6 
m(^V(To)) [top right] by decomposing the indices (a v j,a w j,a Wy 2 , QV,3> i, a^g) into 
pairs. The contribution to the OPE coefficient (!Po)(Eo) = (Co) A f A (y, x) corresponding 
to a perfect matching cr is given according to eq.(4.61) by explicit expressions involving 
the propagator A [bottom]. 

We now want to extend this representation to more complicated trees T. As a 
warm up, let us first consider trees T i with only one internal vertex u besides the root, 
I<r = {u}, such as the tree displayed in fig.l. As mentioned earlier in (3.51), trees of this 

4 The r.h.s. of (4.60) is also called the Hafnian of the matrix X, see [19]. 
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type correspond to a product of two OPE coefficients, (Co)j^* } cch( (Co^j e . Using 
the representation (4.60), we can express this product in terms of two weighted perfect 
matchings: 

z <W'> = z ( x n s -) t x n *») < 4 - 63 ) 

[Au\=D [A U ]=D o- 1 ei>t(W’ 1 1 )) 7 r i <EO 'i a 2 e'm( r V(T^)) ^cr 2 


Here we write T° for the tree which is obtained from T \ by deleting all vertices and 
edges above the internal vertex u 6 I %, and T\ for the tree which results from T\ by 
deleting all vertices and edges below u e !■%, see hg.3. 



Figure 3: The decomposition of a tree T\ into subtrees T“, T 1 ’ without internal vertices. 


Equation (4.63) can be simplified in various ways. Firstly, we note that the internal 
vertex u e 1% appears in both matchings, which we can highlight by writing the above 
equation as follows: 


£ (*Wi) 

[A„]=D 




2 




IT I e®( ( r V(7’“)) [(v,i),(H’,/)]6o-i U tr 2 
o- 2 em( r V(T b l )) v ’ w * u 


z n 


^(v,i),(u,k) 


[A„]=£> [(v,i),(»,*)]6cri 
[(u,k),(w,j)]ecT 2 


• z 


(u,k),(w,j) 


(4.64) 


The product on the very right, which contains all matchings involving the internal 
w-vertex, can then be written as 


X : 

\<*u,k\= d 


- l (v,i)(u,k ) ' u,k)(w,j ) 


X 

I \<*u,k\=d 


d x v i x v Xu) U,k c\ a u,k oPwJ A / \ _ FTpd— |Q'v,ii 

@X U @X„ %U X w ) — 1 v..—»r. A-Iv ; 


O'u.k'. 


y d x l J (x v -x u f u - k d“- (x u -xk) w -j _ 

a uk \ a„, ;! ~ l x v ->x u Hv,i)(w,j) 


V\«u,k\= d 


x v ^x u W ^ R 

w = <R 


(4.65) 


where T d is the Taylor expansion operator 


T Uyf(x) 


\Z M=d ( ^fd;f(y) for d > 0 
10 for d < 0. 


(4.66) 


We can further simplify eq.(4.64) by expressing the summation over the matchings 
ctj, o "2 in terms of matchings <x 6 S0f(X U R(T \)). This is achieved by merging the two 
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Figure 4: Given two perfect matchings cr\ € 2 € '3R( t V(T^)) [top] we obtain a 

matching cr e 9 ft (_£ U K{T 1)) [bottom] by merging the vertices corresponding to the 
internal vertex u e I%(T 1) as indicated by the dashed red lines. 


matchings at the ( u , k) vertices, as shown in fig. 4. Note however that this mapping 
{cr 1 , cr 2 ) —* cr is not one to one: Exchanging two vertices {u, k) and {u,k') yields the 
same matching cr. For a given cr, we therefore pick up a symmetry factor |/(cr)|!, where 
(recall that de(w) denotes the descendants of u) 

/(cr) := {[(v, i)(w, /)] 6 cr : v e d e(u),w £ de(w)} (4.67) 


is the set of merged edges, i.e. those adjacent to a w-vertex in the original matchings 
cr l ,cr 2 . The matching procedure thus leads to the formula 


{ qpkl/r |Of v ,j|yi 

_ n 


if n € /(cr), v e ch(it) 
if n /(cr) 


(4.68) 

where we summarised the possible assignments of the Taylor expansion degrees to the 
merged lines in the definition 


D{cr) - {d - {d n ) neI (a-) : d n 6 N, ^ {d n + 1) - D}. (4.69) 

7T e/(cr) 

We can generalise this strategy to the expression (Pq)(T) for more complicated trees T. 
For this purpose, let us first define the sets 

/«(cr) := {[(v, i)(w,.[)] 6 cr : v 6 d e(u),w £ de(w)} (4.70) 
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for any u £ !<%, which contain all edges which are merged by connecting two n-vertices. 
Further, define 

D(cr) = [d = (dXilr) : < € N, J] « + D = Ad , (4.71) 

neI„(o-) 

which is the set of all assignments of the Taylor expansion degrees to the merged edges. 
We then have the compact formula: 

Lemma 2. LetT eT(x\A). Then 

n i (Po)(T) = e n I4(c0|! zn M n {d) (4.72) 

ueI K (T) [A„]=D„ o-em(£uK(T)) ueI K (T) de£>(<r) necr 


where the matrix M(cf) is given by (recall that by an(v) we denote ancestors ofv) 


M, 


n=[(v,i)(w. 




n 


1 X V ^>X U 


n 


\a w j\ 


(4.73) 


wean(v)\an(w) 


i , ean(w)\an(v) 


The product over the vertices u, s in eq.( 4.73) is ordered from leaf to root, i.e. every vertex 
is to the left of its ancestors. 

Proof The proof works by induction in the number of internal vertices \I<r(T)\. In 
the simple examples above we have already dealt with the cases \I<r(Tq)\ = 0 and 
\I'n(T\ )| = 1, so the induction start has already been taken care of. The induction step 
works as follows: Assuming that lemma 2 holds for all trees T' 6 T(x\ A) with up to n 
internal vertices, \I<r(T')\ < n, we have to show that the lemma also holds for trees T 
with n + 1 internal vertices, i.e. for \I^(T)\ = n + 1. 

The idea of the proof is analogous to the simple example with one internal vertex 
discussed above: Fix any internal vertex u 6 IniT) and denote by T a the tree obtained 
from T by deleting all vertices and edges above the vertex u, and by T b the tree obtained 
from T by deleting all vertices and edges below u. Since both T a and T h have at most n 
internal vertices, we can use the induction hypothesis in order to express fPftfr) as a 
product of the form 

n e (*W) = e e n |/,M! E I"K Ma) 

seI K m [A S ]=D S [A„]=D„ 0 - <! 6®((Xu^(7’“)) seI K (T°) J e2Xo - , n a eo- a 

“ K a> (4 74A 

x z n |/?(crfo)i! z n M ^ b) 

o- b eW(£unTp) S eI K (Tp J b eV(cr b ) n b^cr b 


From here on we can essentially repeat the discussion following eq.(4.63): We dis¬ 
tinguish matchings in cr a and cr b containing the vertex u, and those that do not. For 
the former, we obtain products of the form M (vMu ^ ) (d a )M (ll ^ w j ) (db), which can be 
simplified using eq.(4.65): 

^ ^ M{Yj),(u.k){d a ) • {w jdd /,) — 4/f \ij)(\\!jj(d) (4.75) 

\‘Xu,k\=d% 

Expressing the matchings <r a , cr b in terms of matchings cr 6 9k(X U “R(T)) by merging 
the n-vertices as before (see fig. 4 and the corresponding discussion), we pick up a factor 
|/„(cr)|! and thereby arrive at the representation (4.72) as claimed. □ 


21 


4.1.1 Proof of the bound (3.53) for r = 0: 


Lemma 2 provides a compact expression for the objects of interest in our proof of 
theorem 2. Next we would like to derive an upper bound for the r.h.s. of eq.(4.72). This 
is achieved with the help of the following lemma: 

Lemma 3. Let M(d) be the matrix defined in eq.( 4.73), let b E S(T) be the branch ofT 
fixed in theorem 2, let n = [(v, i)(w, /')] £ cr and define the shorthand 


|(1 + s) foru E b(T) 
|l/e 2 foruib{T). 


(4.76) 


where s E (0,2 ~ T 3 ) with Tv as defined in theorem 2. For any 6 > 0 one has the bounds 


(K,l+Kjl+<5)! n etf"- 4 h( “' nf)r )v 5 +1 

__ utl’n _ 

\M (rf) I < min L«--" c vl) 1+|av ' il (£ 2 min p„-x w |) 1+|ar,, U m s (s 2 \x e -Xf\) s 

\ 1 V 1 7 T\ U J\ — tie sb(v) MEsb(w) 


for v, w ± % 


max„ ech (K) \x u -xrP w L +1 


(4.77) 


min„ esWv ) \x u -xd av ^ +l uePj 


n uerJK-d] 


rCh (u)C\P, 


')•& farw = % 


where we use the shorthand P n := (an(v) \ an(w)) U (an(w) \ an(v)) and where e is the 
vertex closest to the root in the set an(v) \ an(vv) (ancestors ofv which are not an ancestor 
ofw), and similarly for f with the roles ofv and w exchanged. 


The straightforward but tedious proof of this lemma can be found in appendix A.2. 
Using lemma 2 we can bound the l.h.s. of (3.53) for r = 0. 


11 Z (Po)(T) 

uelftiT) [A u ]=D li 


<- £ n |/„(cr)|! zn \mM\ 

o-s»i(£u«) ueI K JeDicr) ne<T 


(4.78) 


We would like to bound the product of matrix entries on the r.h.s. of this inequality with 
the help of lemma 3. For this purpose, we first note that the product of combinatoric 
factors can be simplified using £ ffe /„(o-)(d" + 1) = D u and using 

' |(k>-,,| + \a w j\ + 6)\< (^[A v ])!. (4.79) 

necr ve£ 


A rather non-trivial point concerns the factors m s \x e - xf 6 in the bound (4.77). How 
many of these factors do we obtain in the product over n E cr? Note that, on account 
of the ^-functions, our bound (4.77) on the matrix elements M n vanishes if d £ !D(cr) 
contains two elements d e n < d{ such that e is closer to the root of T than /. Pick a 
vertex u E I(T). If we have A (u) > 0 at that vertex, then there has to be at least one 
pair [(v, i)(w , /)] = n E cr such that v, w E de(w) for the product of matrix elements not 
to vanish, since otherwise we would have a pair d“ < d c ^ (u) . From lemma 3 we know 
that in this case, since clearly v, w ± *R, we have the freedom to generate an additional 
power of 1 / (m ■ min (V ,/ ech(u) \x e - x e > I)* 5 ". Repeating this argument at every internal vertex 
of T, we arrive at the bound 


n i M ^)i < 

necr 


max \x u 

ue ch(«) 


^r* ] (£[A v ])! g~ 2 I,veV\HT)[ A v](l .). gfilweb(T)l A wi 
veX 


n min \x u - x v | [Av] Y\(m ■ min |v, ~ Xj\) 0<A " ) ' 6 " 

v€X «esb(v) „ 6 j i,jech(u) 


I IT ■ 

ieI K 

(4.80) 
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Substituting the bound (4.80) into (4.78) and using also the estimate 


[~[ I4i( cr )|! < n! |J *' 

ueI K 

where n := XvsXu k n v ^ with n v as dehned in (4.59), as well as 

|1(XU^)| = (n- 1)!! < n! 


(4.81) 


(4.82) 


and 

\D{(t)\ < ]~[(D„ + l) n < Y] (D„ + l) n < g _ Zve'VXiUv] n (D ll + lf T , (4.83) 

uel ueb(T) ueb(T) 


to bound the summations over cr e S H1(X U K) and over cf 6 'DUr) in (4.78), we finally 
arrive at a bound for the quantities of interest: 


n z (r " m 

ueI K {T) [A„]=D„ 


< 


max \x u 

uecWR) 




[A«] 


! ] min \x u - v,.| 1 ' 41,1 f| (m • min \xi - Xj\y HA " } ' 6 " 

ve£ «Esb(v) ug j /Jech(«) 


n« 

ueI>R 


x n! |J * l+1 • (Z[ A >’])! £- 3 2v 6 w) [ Av ] p|(Z) vv + 1)^(1 + s) lAw] 

(4.84) 


ve£ 


web 


This inequality indeed implies the bound (3.53) for the free field (r = 0) if we choose 
the constant K such that K > n! |J,?l+1 • (2 v6jC [A v ,])!. 


4.1.2 Proof of the convergence relation (3.55) for r - 0: 

To complete the induction start, it remains to be shown that the convergence property 
(3.55) holds for the free theory, i.e. we need to show that (suppressing for the moment 
the dependence on the coordinates x t ) 

oo 

1™( R o)a 1 ...AmAm + i..a n = (Co)ai...a n - Z (Co)a i ...a m (Co)ca m+ i...a n = 0 ( 4 - 85 ) 

[C]=o 

on the domain £ < 1 defined by (1.11). In terms of our tree notation, we can write the 
associativity remainder as 

= (Pom) - X E (!P " )(7 '' ) • <4 ' 86) 

d<D [C}=d 

where T 0 e T ((A u ... ,A N ,B);(x u ..., x N )) and T x e T((A U ...,A N ,B, C); (x u ..., x N )) 
are the trees shown in figure 1. Using the bound (4.84) for the r.h.s. of this equation, 
one can verify that the sum over d is absolutely convergent on the domain £ < 1 in the 
limit D —> oo [see the discussion following eq.(3.56)]. 


23 




Thus, it remains to show that the limit in eq.(4.85) is indeed zero 5 . To see this, we 
recall equation (4.68), which we can write in the limit D —> o o and for £ < 1 as (using 
the Leibniz rule to pull Taylor expansions out of the product) 


oo 

Z < c »>* ...Am (^Q^CAm+i—An 
[C ]=0 


oo 

= z n -z 


o-€Sm(XU«(7’ 1 ))ff6o-\7(o-) d=0 


^.n ^ 

n'el(cr) 


(4.87) 


Here T'^ XM )-t( XM Vu) is the multivariate Taylor operator, 


T ; 


(x l ,...,XM)->(XM,~,XM)f( XU ■ ■ ■ ’ X M) ~ X 

|vi|+...+|vmI=^ 


M 

n 


(Xj - x M y 

v,-! 






(4.88) 

Using the fact that the Taylor series is convergent on the mentioned domain and 
recalling our explicit formula (4.60) for the zeroth order OPE coefficients, we therefore 
arrive at the relation 


OO 

X (Co )£, ...A m ( C o)cA M+l ...A N 

[C]=0 


crem'VtTo)) neo- 


(4.89) 


which establishes equation (3.55) for the free field and thereby concludes the induction 
start. 


4.2 Induction step: Higher perturbation orders 

Assuming that theorem 2 holds up to perturbation order r, we now want to show that 
it also holds at order r + 1. Our main tool to achieve this task is the recursion formula 
for the OPE coefficients, eq.(1.12), which in turn implies a corresponding recursion 
formula for the expressions ( P,)(T ). 

4.2.1 Proof of the bound (3.53) at order r + 1: 

When expanded in g, our recursion formula 6 (1.12) reads at order g r+1 : 

0 Cr + i)i.. AN (xi, ...,x N ) = —^ J d 4 y [(C,.)£ Al Ajv (y, x u ..., x N ) 

N 

X (Cn ’ Xi) ’ ' ' ' ’ 

i=l [C]<[A,-] ri+r 2 =r 

- X X ( C n)A l ...A N ( X ^ ■ • • >X N ) C Cr 2 )l c (y, Tv)] , 

[C]<[B] ri+r 2 =r 

(4.90) 

5 This fact has been shown previously, in [12] for the case r — 0,N = 3. 

6 Our choice of “renormalisation scheme" enters at this stage: The particular form of the recursion 
formula given here was derived for the so called BPHZ renormalisation conditions. See section 5 for a 
discussion of renormalisation ambiguities. 


24 






where the index 2 corresponds to the interaction operator of our model, i.e. 0& := (f 4 /4!. 
Formula (4.90) allows us to write the l.h.s. of (3.53) at order r + 1 in terms of r-th order 
quantities via 


n z^ r )- r+1 n x x 

i£.I<R [Aj]=Dj [Ai\=Di Yu U £i r u =r +1 veX 

= ^n Z Z f [ (P ' )(7V) _ Z E ~ E 

ieI K [A,]=A veJ ^ wech(v) [A„]<[A W ] [A„]<[A V ] 

(4.91) 

where the trees T v , T vAu , T Au V are obtained form T e T(x\ A) as follows (see fig.5): 

• T v 6 T ((-?,)’); (A, 2)) is obtained from T by connecting an additional leaf with 
weight (2,y) to the vertex v. 

• T vAu g T{(x,y)\ (A, 2, A u )) is obtained by connecting a leaf with weight (2,y) 
to the parent edge of v, splitting this edge into two halves. The new vertex u 
adjacent to these two halves receives the weight (A„, x v ). 

• Ta u , v 6 7”((f,y);(A,2,A„)) is obtained by connecting a leaf with weight (2,y) 
to the parent edge of v, splitting this edge into two halves (if v = R, then we 
add a parent edge to v and connect the leaf to this new root). The new vertex 
u adjacent to these two halves receives the weight (A v , x v ), and we change the 
weight of the vertex v to (A„, x v ). 



Figure 5: The trees T v , T vAu , T,\ ij V are obtained from the tree T by adding an external edge. 

Our plan is now to combine the formula (4.91) with the inductive bound (3.53), which 
holds up to order r by assumption, in order to verify the bound (3.53) at order r + 1. The 
terms under the integral in eq.(4.91) can be estimated with the help of the following 
bounds: 
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Lemma 4. Denote by B r (T) the r.h.s. of (3.53). Then 


n z (p ' m 

izlft [Aj\=Dj 

(Y\web(D w + 1 ) 

S\ 


< 


B r (T)K 

min \x w - v|‘ 


sup 


\y - x v 


M'Ech(v) 
4-8 r+1 

sup 


min \Xj — Xj\ 

<'Jech(v) 

min \y - a,-| ’ max \x { - x , 

y (Ech(v) i6ch(v) 


r,l 


2 <5„ 


n 

/ech(v) 


f min |a ; - - jc,- \ 

;'€ch(v)\{i) 

— ? * 


ly - */l 


[A/] 


sup 


/ 


ly - 


max |A', - Xv 

Vjech(v)\{Z) 


[A v ] 


(4.92) 


Z n Z <^xr vA ) 

[A„]=a [A/i—ii/ 


z n z < p ' xT ^> 

[A„]=a (e/r [a,]=a 


< 


fi r (r) K sup(l, l/m|A v - y|) 




jfl^([A v ]+4-d)6 u — y|4+0([A v ]+4— d)8 u 

( min \xi - Xv\ 4 5pa(v) 

x sup 


min | A y - a, 

i’esb(v) 




l V - Ay 


[A„]-a 


i€sb(v) 


ty-* v 


i 


n«€/)(^v)> 4-1) 


4-8 r 


< 


B r (T) K sup(l, l/mlA v - yl)* 5 " 

YflO(d+4—[A v ])6 u \jCv — y\4+6(d+4— [A V ])<5 M 


ly - Xv 


X (v,d) (4.93) 

\[A„]-d 


max | a v - a, | 

y ('ech(v) 


( min \Xj - Xj\ \ 


x sup 


»t/Ech(v) 


l y - Xy\ 


1 


A 


Wweb^Pw "4 1 ) 


4-8 r 


X(v, d) 


(4.94) 


where K > 0 is a constant that depends neither on the integers D\ nor on s or the 6 V , and 
where we use the shorthand 


X(v, d) 


j(l + e) 8 " 1 d ■ (d + l) 8 ' +1 (®r+4) z y v e b 

jg~ r+lrf ifvib. 


(4.95) 


Proof. The lemma follows by straightforward computation from the inductive bound (3.53). 

□ 


We now substitute these bounds under the integral in the recursion formula (4.91). It 
is, however, not possible to estimate the resulting individual terms directly, because the 
integral over each individual term contains divergences in the regions where y w a; 
(UV) or where |y| » sup ; |a,| (IR). As mentioned in our overview of the proof at the 
beginning of this section, we therefore have to take a little more care and take into 
account cancellations between these divergent terms for each of those dangerous 
regions. In order to study these cancellations of singularities at short- and large 
distances, we define the following partition of R 4 : 

Definition 5 (Integration regions). Let v £ I(T) be an internal vertex of the tree T and 
let b £ !B(T) be the branch mentioned in theorem 2. Then 

(UV-regions) 


26 

































ny: 


{y e R 4 : \Xj - y\ ■ (1 + e) 2 ' 8 ' + ' < min 7 - €sb(0 \x t - jc ; |} ifi G ch(v) n b(T) 
{y G R 4 : | Xj - y\ • g - 2'8 r+1 < min 7€sb(0 \x t - *,]} ifi G ch(v) \ b(T ) 


(4.96) 


(IR-region) 

{y G R 4 : \x v - y\ > max \x v - xj\ ■ (1 + e) 2 ' 8 ’ +1 


y'Ech(v) 


}\u ; ny if pa(v) G b(T) 


n v IR : = 



x v -y\> max |jc, 

y'sch(v) 


; v - Xj\ • £ 2 ' 8r+1 ) \ U,-QV if pa(v) ^ &(r) 


(4.97) 


(Intermediate region) 

n.' IM ;= R 4 \ (U,Q) u a v IR ) 


(4.98) 


Remark: Note that for any v G 1(7 ) one has Q.'j M U if R U,g ch ( V ) If = R 4 and that 
these sets are disjoint, in particular If fl D 1 . = 0 if / y j. Note further that the UV- and 
IR-regions get smaller as we increase the perturbation order, which will be needed later 
in order to obtain sufficiently strong bounds within those regions [more precisely, this 
fact is going to be crucial for the estimates (4.106) and (4.112)]. 

We now derive a bound on the r.h.s. of (4.91) by decomposing the y-integral into 
integrals over the regions defined above. We will see that, indeed, the resulting bounds 
for the contributions from each of those regions are consistent with (3.53) at order 


r+ 1. 


The intermediate distance region In this region the integration variable y of 

eq.(4.91) is neither very close to, nor very far from the points (jq)iEch(v)- Hence, we will 
encounter neither UV- nor IR-divergences, and we can simply insert the bounds from 
lemma 4 in order to estimate the integrand, without taking into account any further 
cancellations. 

Let us fix an internal vertex v G I(T). By definition, we then have for any e G ch(v) 



(4.99) 


Furthermore, we have the inequality 



k/-*vl 

e 2-8 r+1 


7M 


(4.100) 


' max | Xi - x v 


< (2 n) 2 n 


i€ch(v) 


e 4 ' 8r+1 min \x t - x\ 

V y'Ech(v) 
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where as before n := Xve£u« n v is the “total number of fields" associated to the external 
vertices of the tree T. Combining these inequalities with lemma 4 and choosing 8 
sufficiently small such that 8 + 8 V < 1, we obtain for the first term under the integral in 
(4.91) the bound 


Jinx 

0 V ieI K [A,]=Z 


(Pr)(T v ) 


d 4 y 


< B r (T) K 


/rwAv + m 

/ ^ ^ £^2ee(ch(v)Uv)nfo[^e] \ 

' max | Xi - x v \ ' 

(Ech(v) 

\ s l+6 ) 

1 ^^je€(ch(v)Uv)\fe[^e] I 

min | xi - xj\ 

y i*je ch(v) / 


(4.101) 


where constants (i.e. factors depending neither on the weights D, nor on e) were 
absorbed into K. The last factor on the r.h.s. can be absorbed into the expression B r (T ) 
by adjusting the parameter 8 V —> 8 V + 8 e (0,1). To see that the resulting bound is 
smaller than the r.h.s. of (3.53) at order r + 1, we note that the inductive bound (3.53) 
grows like 


B r+ i(T) = B r (T) K 


(1 + e)^ H 


e,[A w ] 


^2ve'V\h[-^v 


\7-8 r+1 

• ^ 

web / 


(4.102) 


as we increase the perturbation order r, where K is some constant that depends neither 
on the Dj nor on s of the 8 V . Since the remaining terms on the r.h.s. of (4.101) are indeed 
smaller than the r.h.s. of (4.102) (choosing 8 < 1/4 and assuming that XveV^IA] > 0), 
we conclude that this contribution to the recursion formula (4.91) is consistent with 
the claimed bound (3.53). 

Similarly, using lemma 4 as well as the estimates (4.99) and (4.100), we obtain for 
any w e ch(v) the following bound on the second term under the integral in (4.91): 


fi x nz (P r )(T wAu )\d 4 y<B r (T)K 

rt vn r/t.n-n. 




X 


A U \<D W ieI<R [Ai]—D, 


( max(l /m, \x { - v,. )) 

3<5„ 

i'ech(v) 

X | 

min | Xi - X;\ 

V iVisch(v) , 



FLe/XA + 1) 


r.2<5„ + 1 


4-8 r 


\{D W + l) 8r+1 (®r+ 4) + 1 (l + 6 .)2- 8r+1 (A,+5 v) 
I (D w +1) e -2( D w+'5v ) 8 r+1 


if w 6 b 
if w i b 

(4.103) 


The factor with exponent 8 U can again be absorbed into B,(T) by choosing 8 U sufficiently 
small and increasing the value of 8 V slightly. One checks, using also the inequality 
(. D w + 1) < s~ Dw for the case w £b, that the bound (4.103) is indeed smaller than (4.102), 
and it is therefore consistent with our hypothesis (3.53). For the third term on the r.h.s. 
of (4.91) we can proceed in essentially the same manner as for the second one and we 
find that also the integral over | X[a„]<z>, Xhei K X[a,]=d,(^-)(^a,„v)I satisfies the bound 
(4.103). 

Thus, we have found that the contributions from each term under the integral are 
smaller than the claimed bound (3.53). In order to bound the total contribution from 
this integration region, it remains to estimate the number of terms appearing under the 
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integral. For a given v, the integrand in (4.91) contains | ch(v)| + 2 < n + 2 terms. The 
sum over all vertices v contains \I(T)\ < n terms. Both of these factors can be absorbed 
into the constant K in our bound. 

To summarise, we have verified that the contribution to the r.h.s. of recursion formula 
(4.91) from the intermediate integration region is smaller than the claimed bound (3.53). 


The UV-regions Q'.: Here the integration variable y is close to one of the points 
Xi, so we have to take into account cancellations between different terms under the 
integral in the recursion formula (4.91). In order to achieve this, we not only make use 
of the inductive bound (3.53) here, but we also apply the induction hypothesis (3.55) 
stated in theorem 2 in order to organise the short distance cancellations. 

Fix a v 6 I(T ) and a we ch(v) and consider now y G QJJ,. To bound the integral 
over the expressions (P r )(TA u , v ) and (P r )(T Au )with i e ch(v) \ in’}, we can proceed as 
above and arrive at the same bounds as in the intermediate region y G fl) vr For the two 
remaining terms under the integral, our second induction hypothesis, eq.(3.55), implies 

<!Pr)C Tv)” 2 ( P r)(T wAu ) = (P r )( T wAu ). (4.104) 

[A u \<n w [A u \>n„ 


To bound the r.h.s. of this equation, we now use lemma 4, distinguishing the cases 
w G b(T) and w i b(T ) in the process. Making use of the inequality 


|y - x w \ < |(1 + e) 2 ' 8 ' + ‘ if w e b(T) 
min \x w - Xi\ ~ | s-* r+> if w t b{T) 

z'esb(w) V 


for y G Q. 


V 

w ’ 


(4.105) 


we obtain the bound 

1 1 z n |dv 


[A U ]>D W ieI<R [Aj]=Dj 


J r* o o 

d4 >' z 

a::. 


<5,, 


n^([A/] + i)\ 48r+1 


x 


B r (T) K max(l, l/m|y - x w \) 

,X W - y\3+OQ-d)6u m O(3-d)6 0 m j n _ x .| 

^ z’esb(w) 

[(1 + £) 8 ' +1 W»- d+l Xd + D w + 2) sr+1( ^+ 4 ) for w G b(T) 


-8 r+1 (D w -d+ 1) 


< B,fT) K 


n^([A«] +1) 


4-8 r 


sup(l, 


for w i b(T) 

\2 5 U 


m ■ min \x w - x { \ 

i'esb(vk) 


X 


1(1 + sf +l (D "- +2) £-2-8 r+1 (33r+4)-2 (gr+l^ + 4 ))| fo r w 6 b{T) 


-8 r+1 (A,+l)-l 


for w i b(T ) 


(4.106) 


Here we used the inequality 


oo , \(D +2)8 r " 

£« + + S + 4 »> 


d =0 


(4.107) 
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as well as the elementary estimate 


z 


£ 


\ r+l d 


< 


1 1 

- < - 

1 - £ £ 


(4.108) 


to bound the infinite sums and the inequality (choosing 8 U < 1/2) 



max(l, l/m|y - x w \) Su d 4 y 2 / 1 

... -775—777-:-:-: < (27 T) SUO 1, -:-;-; 

y\i+vv-a)d„ m v(i-a)o II mm _ Xj \ \ m . mln | Xh , - jq| 

i€sb(w) i'esb(w) 


(4.109) 

to bound the y-integral. Choosing <5„ small enough such that 8 V + 2b„ < 1, we can 
absorb the factor with exponent 2 8 U into the bound B r (T) via a redefinition of 5 V . The 
factorial (8 r+1 (®7' + 4))! can be absorbed into the constant K. As the remaining terms 
on the r.h.s. of (4.106) are smaller than (4.102), we conclude that also this contribution 
to (4.91) is consistent with our inductive bound (3.53). 

To summarise, we have verified that contributions from the integral over the short 
distance regions Qj’, to the r.h.s. of (4.91) are smaller than the claimed bound (3.53). 


The IR-region Q' IR : Here the integration variable y is far away form the points ;q, 
and we again have to take into account cancellations between different terms under 
the integral in order to bound this contribution to the recursion formula (4.91). 

Fix a vertex v £ I(T). For the second term on the r.h.s. of (4.91) we can proceed 
essentially as in the case of y 6 Cl v IM before. The only difference here is that instead of 
(4.100) we use the inequality 


f 

IR 


d 4 y 


< 


(2tt) 2 


m s "|y - x w \ A+6 “ (m min,- esb ( w) \x w - V/l) 15 " 


(4.110) 


to bound the integral over y. This factor can be absorbed into a redefinition of S v as 
explained previously below (4.109). 

In order to find useful bounds on the remaining terms, we again have to make use 
of our second induction hypothesis, eq.(3.55), which implies that 

C Pr)(T v )~ (P r )(T Au , V ) = J]( P r)(T Au J (4.111) 

[A„]<A, [A„]>A, 


for y £ Q.' ir . Lemma 4 then implies for the r.h.s. 


X. I z n £ 


d 4 y 


IR [A„]>Ov ieI K [A;]=A 


J r> oo 

d *y Z 

Vr T~Tx 


B,(T) K max(l, \/m\y - x v \)° u + 1) 

\x v - y\ A+6 «m Su 


4-8 r 


X 


< 


[ (1 + £f +l (2D '- d \d + l) 8r+1( ® r+ 4) for V 6 b(T) 
|^ g -8 ,+1 (2 D v -d) for V i b(T ) 

r.h.s. of (4.106) 


(in ■ maxy 6ch ( v) \x v - x fi 6 - 


(4.112) 
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Here we used the same estimates as in the short-distance case to bound the sum over d, 
and we used 


f 

Jyei 


d 4 


>’ 


(2 nf 


M (4.113) 

kv - >f +,i max ;ech(v) \x v - xj\ 6 

to bound the ^'-integral. Choosing S u small enough, we can absorb the first factor on 
the r.h.s. of (4.112) into a redefinition of 5 V . The remaining terms in the bound (4.112) 
are then smaller than (4.102), and we conclude that also this contribution to (4.91) is 
consistent with the claimed bound (3.53). 


Combining our bounds for the intermediate-, UV- and IR-regions, we conclude that 
the r.h.s of (4.91) satisfies a bound that is smaller than our hypothesis (3.53) at order 
r + 1. 


4.2.2 Proof of the convergence relation (3.55) at order r + 1: 

The last step in the induction is to show, assuming that theorem 2 holds up to pertur¬ 
bation order r, that the second statement of the theorem, eq.(3.55), holds also at order 
r + 1. For this purpose, we write down the recursion relation for the remainder, i.e 
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where T 0 e T ((A u ..., A N , B); (x u ..., x N )) and T x 6 T((A U ...,A N ,B, C); (x u ..., x N )) 
are the trees depicted in fig.l. In order to show that this expression vanishes under the 
assumption £ < 1, we would like to exchange the order of the integral and the limit. By 
the dominated convergence theorem, this is allowed under the following conditions: 

1. For all D 6 N the integrand is bounded by some integrable function B(y). 

2. The limit D —* oo of the integrand converges pointwise almost everywhere. 

The first condition is easily checked with the help of the bounds derived in the previous 
section combined with the inequality (3.57) to bound the sum over [C]. For the 
bounding function B(y ) we can choose for example 
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for some 6 £ (0,1) and for s £ (0,2 _x>r_4 ' _3 ], where X> r = £,[A ; ] + [5], To show 
that the integrand converges pointwise to a limit as D —> oo, we make the following 
observations: Using our induction hypothesis (3.55) at order r, it immediately follows 
that the last two lines of (4.114) vanish in the limit under the assumption £ < 1. To treat 
the remaining terms under the integral, we have to take a little more care: Consider 
first the region 

:= {y 6 R 4 : \x M -y|(l + e) 2 8+1 < min \x M - x\ , |y - x\ > 0} (4.116) 

M< j<N 

for some small s > 0. In that case, the first two terms under the integral in (4.114) 
cancel in the limit D —» oo by our hypothesis (3.55), and the remaining terms under the 
integral are of the form 
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lim B r+ \(Ti) = 0. 


The first equality follows simply from eq.(3.55) at order r, and the estimate in the third 
line follows analogously to our discussion of the short distance region in section 4.2.1 
[see (4.106)]. Thus, we find that for yefij the integrand converges to 0 as D —> oo. 

In the region 

f2 2 := {y 6 R 4 : \x M -y| > (1 + e) 2 ' 8 + ' max \x M - x;\ , | y - x,| > 0} (4.118) 

1 <j<M 

we simply exchange the role of the second and third term on the r.h.s. of (4.114) and 
otherwise proceed in a similar manner, using estimates from the previous discussion of 
the large distance region Qjr [see (4.112)]. We find that the integrand also vanishes in 
this region. Note, using the assumption £ < 1 and choosing e sufficiently small, that the 
two regions and Q 2 cover all of R 4 apart from the zero measure set [y = jq, i < N}. 
Thus, we conclude that the integrand converges pointwise to 0 almost everywhere. 

To summarise, we have verified that we are allowed to exchange the order of the 
integral and the limit in (4.114). Since the integrand vanishes in the limit, the same 
is true for the integral, which establishes the second statement of theorem 2 at order 
r + 1, thereby closing the induction and finishing the proof of theorem 2. □ 


5 Massless fields 

The associativity proof for the OPE presented in section 4 was restricted to the case of 
massive fields, m 2 > 0. In fact, the main ingredient in our construction, the recursion 
formula (1.12), only holds for massive fields as stated. In the naive massless limit, the 
right side of the recursion formula becomes ill defined already at first order in g. This 
feature, however, does not indicate a fundamental problem with our approach, but 
is basically due to the fact that our definition of the composite operators (implicit 
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in our recursion formula) is unsuitable for m 2 = 0 . To get around this, we will first 
apply a field redefinition (for nf > 0 ) as introduced in definition 2 of section 2. A field 
redefinition changes the OPE coefficients as in eq.(2.27). Consequently, these will also 
satisfy an appropriately modified version of the recursion formula (1.12). It turns out 
that a field redefinition (depending on an arbitrary “scale” L > 0 ) can be found such 
that the modified recursion relations possess a well-defined limit nr —* 0 . At this stage, 
the same procedure as in the massive case can then be applied to prove the associativity 
property claimed in theorem 1 also for massless fields. 

5.1 Recursion formula for massless fields 

Consider a field redefinition in the sense of definition 2, which is written in terms of a 
mixing matrix Zf e C[[y|] as 

O a = Y j Z B a ° b , (5.119) 

B 

where 0 A are the redefined fields. The matrix Z B has to be invertible in the sense of 
formal power series and it has to be “upper triangular” in the sense that Zf = 0 for 
all [. B ] > [A]. The corresponding transformation for the OPE coefficients is given by 
[compare (2.27)] 


yrB _ \ 1 \ 1 yC i yCN/y—l\B s~>Cq 

^Ai..A n ~ /_j " ' /_j A4j " ' AVv'A /Co ^Ci...Cjv ’ 
Co Cn 


(5.120) 


where we note that all summations are finite because Z is upper triangular. Combining 
eqs.(1.12) and (5.120), we immediately see that the redefined OPE coefficients now 
satisfy the recursion formula (suppressing spacetime arguments) 
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(5.121) 


where the objects T^ are defined as elements of the matrix T, 

T := -Z -1 d g Z . (5.122) 

We would like to make a specific choice of the mixing matrix Z in order to cancel the 
contribution to the integral (5.121) coming from large ]y| (infra-red region). For that 
purpose, we define 



|0 
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for [A] > [ B] 
for [A] < [5] 


(5.123) 


for some L > 0. (Note that T depends on g .) The solution to eq.(5.122) can then formally 
be written as 

Z(g) = V exp - fW, (5.124) 

Jo 

where J J cxp denotes the “path ordered exponential". 

Combining this definition of Z with (5.121) and with the associativity property 
(1.10) and choosing L > max,- |v, - x v |, we can rewrite the recursion formula for the 
new OPE coefficients C B , , as 

A\...An 

d A...A N (* i. ...,x N ) = - f d 4 y ^ (Z-‘)f [c| Al ... Aiv (y, ..., x N ) 

\y-x N \<L ^ 

2 ^A i Xj,^)C^ 1 ...D...A^l-"^w) ^ 5 - 125 ^ 

«=1 [D]<[A,] 

— X? Al Ajv (xi,..., x N )C ED (y, .yv)] • 

|0]<[«l 

Here the idea behind our redefinition (5.123) becomes apparent: We have arrived at a 
modified recursion formula which includes only integrals over a region of finite volume. 
The contributions to the integrals from |y - .y, v l > L have been cancelled precisely by 
the terms coming from the field redefinition (using also the associativity theorem 1). 

We would finally like to tidy up the factors (Z ')(; in front of the OPE coefficients 
in (5.125) by a redefinition of our coupling constant g. In particular, we would like to 
choose this redefinition of g in such a way that the formula (5.125) has a simple and 
well defined limit m 2 —* 0. The following lemma allows us to understand the small 
mass behaviour of the mixing matrix Z: 

Lemma 5. The mixing matrix behaves as 

lim [Zl ■ = 6$ + K ■ 6^ (5.126) 

mr —>0 

for some formal power series K(g). 

Proof We establish this lemma by analysing the small mass behaviour of the OPE 
coefficients appearing in the matrix elements Zf More precisely, we will prove that 

J' (Cr)i. £ (x) = K a ■ [log(L 2 w 2 )] r+I + C>( [log(L 2 m 2 )]' ) for A : 0 A 6 {<p 4 ,(pd 2 p} (5.127) 

\x\>L 

J (C r )Ux) = 0( [log(L 2 m 2 )]' ) for [A] < 4, 0 A t {p 4 , <pd 2 p} (5.128) 

\x\>L 

for some constants K A which depend on the perturbation order, and where K A ± 0 for 
0 A = p 4 . These equations then imply that the rescaled matrix Z;) /Z'~ vanishes in the 
limit m 2 —* 0 unless O a - ip 4 or 0 A - (pd 2 tp, which upon inversion of this matrix leads 
directly to the lemma. 
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To prove these statements, we are going to proceed inductively. Using eq.(4.60), 
one checks (5.127) and (5.128) for the free theory by straightforward computation. For 
the induction step we make use of our original recursion formula (1.12). Using the 
associativity property (1.10), we can rewrite the recursion formula in the useful form 
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where the regions Q, C R 4 are defined as 
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(5.130) 
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for some s > 0. Note that the infinite sums in eq.(5.129) are absolutely convergent by 
our theorem 1. Considering first the case A = £ and focusing on the contributions of 
leading order in log(m 2 ), we are left with 
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Here we used the induction hypotheses, eqs.(5.127) and (5.128), in order to estimate the 
small m behaviour of the coefficients C(!„ and we used the bounds 
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for [A] + [5] - [C] > 4, and 


n\x |=A 

< 

CW=L 


(C r ) AB (x) d 4 v < a 4+[C]-[A]_[B] c( [i og(A 2 m 2 ) j' ) (5.135) 


for [A] + [B] - [C] < 4 in order to estimate the other OPE coefficients appearing 
in (5A29). These bounds can be established inductively: They are easily verified at 
zeroth order using eq.(4.60), and, using the recursion formula in the form (5A29), one 
picks up an additional power of log(L 2 m 2 ) with every iteration. Furthermore, we 
also used the fact that g, ^ = 0(m 2 ) to obtain (5.133), which can also be shown 
inductively using d 2 A{x) = -nrA(x) + S(x). Applying the induction hypothesis (5.127) 
in order to estimate the remaining terms in eq.(5.133), we see that indeed we obtain a 
non-vanishing contribution of the order [log(L 2 m 2 )] r+2 , as claimed. 

The other estimate stated in eqs.(5.127) follows directly from (5.134). Regarding 
(5.128), we note that the zeroth order OPE coefficient vanishes. Using this in 

the recursion formula (5.129), one can verify (5.128) by induction. For the integral over 
the coefficients C 4 a , with [A] < 4 one can even check that the limit nr —» () is finite at 
zeroth order, so (5.128) certainly holds at higher orders by iteration of the recursion 
formula. □ 


Combining lemma 5 with a redefinition of the coupling constant 


(5.136) 


we arrive at the following 

Proposition 1. There exists afield redefinition, eqf 5.119), and a redefinition of the 
coupling constant, eqf 5.136), such that the recursion formula for the redefined OPE 
coefficients has a well defined massless limit. Form 2 = 0 this formula reads 

d A..A N (x i, ...,x N ) = - J d 4 y [c| Al ... AjV 0\ x x ,...,x N ) 
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;=i [D]<[a,] 

[Z>]<[B] 

for any L > max,-1 x, - xfi in the sense of formal power series in g. 

Proof. Using lemma 5 in eq.(5.125), it only remains to show that the contribution from 
the sum over E with Oe = pd 2 p vanishes. This is achieved by induction. Using 
eq.(4.60), which also holds for the new coefficients since (Co) = (Co), and using also 
the fact that (d 2 + m 2 ) A(jc) = b(.r), one verifies that the term in question, i.e. 
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vanishes at zeroth order in the limit m 2 —> 0. To show that this term also vanishes to 
all orders in perturbation theory, we write the corresponding recursion formula in the 
form 7 


<9g[(5.138)] = - J d 4 z (Z" 1 ^ [T ^ v(y) [C B EAl .. AN (z, x u ..., Tv)] 

\z-x n \<L ™ 

N _ 
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(=1 [C]<[A,] 

y ^ r ^-<fd 2 <p(y)\C Al An (X{, . . . , T\) CJ i A '(Z.. Tv)]J 
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(5.139) 


where we defined the operator 


(5.138) =: T^ (y) [C B Al ^ AN (xu -. • ,Tv)] (5.140) 

which acts on products of OPE coefficients by the Leibniz rule. Thus, assuming the 
expression (5.138) vanishes up to perturbation order r, it follows from eq.(5.139) that it 
will also vanish at order r + 1. This closes the induction and proves eq.(5.137). □ 

One may view eq.(5.137) as providing a definition for the OPE coefficients of massless 
</> 4 -theory: We simply define the OPE coefficients of the massless theory to be the 
obvious ones in the free theory [i.e. setting m 2 = 0 in eq.(4.60)], and then define the 
higher orders via eq. (5.137). The OPE coefficients of this massless theory are then 
defined as a formal series in g. 

5.2 OPE associativity for massless fields 

Defining the OPE coefficients of the massless theory via proposition 1 as discussed in 
the previous subsection, the theorem is again that the resulting definition is consistent, 
i.e. does not lead to UV-divergences at any order and satisfies the associativity condition 
at any order in g: 

Theorem 3. The OPE coefficients of massless Euclidean ip A -theory, as defined recursively 
through eqs. (4.60) and (5.137), satisfy the associativity property (1.10) on the domain 
(1.9) to any order in perturbation theory. 

Sketch of proof: With the modified recursion formula (5.137) at hand, we can copy our 
strategy from the massive case in order to prove associativity of the OPE also for 
massless fields. As the differences in the proof are minor, we refrain form repeating the 
lengthy calculations here. Instead, we only point out the main adjustments that have 
to be made. 

7 In the derivation of (5.139) we have exchanged the coefficient for the coefficient . 

This is a non-trivial procedure in the case where C e {(dip) 2 , (ipd 2 ip)}, since in that case these coefficients 
do not actually coincide. However, we note that in (5.139) they multiply vanishing contributions of the 
form Tc(v)[C^ An (x\, ..., Xjv)], so exchanging the order of the indices is indeed justified. 
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Most importantly, one has to adapt the induction hypothesis (3.53) to the massless 
case by replacing factors of 1/m by the length scale L appearing in the modified 
recursion formula. The induction step remains largely the same. Here we can simply 
take the limit m —* 0 in the bound (4.84), which forces us to choose <5„ = 0. The only 
essential difference appears in the estimation of the recursion integral (4.91) over the 
large distance region Q, IR . With the modified recursion formula, this region now has a 
cutoff L. The estimates (4.110) and (4.113) are therefore replaced by 



< (2tt) 2 


miflzesb(w) l* w X t \ 


(5.141) 


for any 6 > 0. Taking into account these adjustments, the proof carries over from the 
massive case without further complications. □ 


6 Conclusions 

In this paper we have shown that the operator product expansion in Euclidean <//j- 
theory satisfies an associativity condition that was originally conjectured in [1], The 
model is therefore the first non-trivial example of a quantum field theory satisfying 
all the axioms of the framework proposed in [1] (see also sec. 2 of the present paper). 
Further, all results derived in that paper which were based on the assumption of 
associativity, i.e. the coherence theorem, the formulation of perturbation theory in 
terms of Hochschild cohomology and the relation to vertex operator algebras, are now 
established within Euclidean ^-theory as a corollary of the associativity theorem. As 
a side result of the present paper, we have also shown how to adapt the recursion 
formula for OPE coefficients, which was originally only derived for massive fields, to 
the massless case. 

The method of proof followed in the present paper can be straightforwardly adapted 
to other self-interacting Euclidean quantum field theory models. Hence, the associa¬ 
tivity condition should also hold for example in the Euclidean Thirring- and the 
Gross-Neveu model. 

Generalisations of our result in various directions would be of interest, e.g. to 
theories with gauge symmetry or to models on curved background manifolds. In 
particular, it may be possible to generalise the finite volume recursion formula (5.137) 
to (Riemannian-) curved manifolds if the scale L is chosen small enough such that one 
can use Riemann normal coordinates to study the y-integral. By far the most exciting 
potential application of our results is that they may help to give a non-perturbative 
definition of quantum field theory in the sense outlined in section 2. 
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A Zeroth order bounds 


Below we derive explicit bounds on zeroth order OPE coefficients which are used to 
verify the inductive bound (3.53) at the induction start r = 0. More specifically, we first 
estimate Taylor expansions of the Euclidean propagator in section A.l and then apply 
the resulting bound in section A.2 in order to verify the estimate claimed in lemma 3 
above. 


A.l Taylor expansions of the propagator 

For free quantum fields, the operator product expansion is closely related to the Taylor 
expansion of the propagator. As we have seen for example in lemma 2, the same holds 
true for the contractions of OPE coefficients Pq(T) considered in this paper. It should 
therefore not come as a surprise that a central ingredient in our derivation of the upper 
bounds on \Pq{T)\ are bounds on Taylor expansions of the propagator. More precisely, 
we make use of the following lemma [recall that by A(v) we denote the Euclidean 
propagator, eq.(4.62)]: 

Lemma 6. For any s 6 (0, ^], any 6 e [0,1], any w £ N 4 and any (d\,...,d r ) e N r , one 
has 
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Proof. Our strategy is to pull the modulus into the summations as follows, 
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To bound the last factor on the r.h.s., we write it as a contour integral: 
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(A. 144) 


Here we use the shorthand u := Vi + ... + v r _i + w, and y is any circle around the origin 
in the complex such that d u A (y + zx r ) is holomorphic on the closed disk bounded by 
this circle. Since the propagator has a pole at the origin, y is restricted to circles with 
radius R < |y|/|v r |. We therefore write 
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where s > 0 is arbitrary. From eq.(A.144) we then obtain the bound 

sup zey \d u A(y + zx r ) | 
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In order to estimate the numerator, we write the propagator explicitly as 


i r 

= - d u 

167T 2 Jo 


d ll A(x) 

Using the inequality [6, eq.(56)] 
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Substituting this estimate in (A.146) and noting that sup_ sy (l/|y + zx r |) = (1 + e)/(e|y|), 
we arrive at the bound 

\ dr / r\ / t \ \ |w|+(5+2 
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Combining this bound with the inequality 
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and choosing s < d- we finally arrive at the claimed bound (A. 142), which finishes the 
proof of the lemma. □ 


A.2 Proof of Lemma 3 

We want to derive a bound on the matrix elements M n defined in eq.(4.73), where 
n - [(v, i)(w, /')] 6 cr for some perfect matching cr 6 U R). Let us first assume that 
v, w ± R. Further, let us write explicitly an(v) \ an(w) = (u\,... ,u a ) and an(w) \ an(v) = 
(si,..., Sb), where we use the convention that w,- is closer to the leaves than u l+ \ , and 
the same for Sj. We can then write equation (4.73) explicitly as 
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Using the bound (A. 142) from lemma 6, we obtain 
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for any £ 6 (0, l/ 8 (a + b)\. Since 8 (a + b) < 8 |J^| < 8 |X| < 8 T 7 - ^ 2 T ’ r+3 , we can 
always choose e 6 (0,1 /2 Sr+3 ], which already establishes lemma 3 for the case where 
{u u ...,u a ,si,...,s b } C\b(T) = 0. 

Thus, assume now that one of the vertices m, is in b(T). In this case, we note that 
also the vertices Uj+i,..., u a belong to b(T) on account of being ancestors of w,. Further, 
also know that none of the vertices (s^,. .., s b ) belong to b, since none of them is an 
ancestor of u { by definition. If u a _! 6 b(T), then it is easy to see that the sum over a Ua l 
simply yields a Kronecker delta 6 d «a d « n \ since by definition all vertices in b have the 
same associated coordinate, i.e. x Ua = x Un , in that case. We can repeat the procedure 
with the line u a -2 if it is in b(T) as well. Renaming summation indices, we can therefore 
reduce (A.152) to a form where only the index u a corresponds to a line in b(T). Thus, 
we see that vertices in b(T) come with factors of (1 + s) instead of 1 /e 2 , which is also 
consistent with the bound (4.77) in lemma 3. 

Next we come to the case w = K. In this case an(vv) = 0, so M„ is simply 


\M n \ = 


x 


< 


< 


£ - Z 

\a H \<<fy -|a v ,i| \<*u a \<d u „ a -\a v ,i\ 

(x v - (x ut - x U2 ) a “ 2 - a “i {x Ua _ t - x Ua ) aua ~ au “- 1 (x Ua - x< R ) awJ ~ av ’ i ~ aua 


a Ul \ 


(Q’m? ®«i)' 


(Q'Ma ^U a - 1)’ (®wj (X\\l OC Ua )! 


kv - Xu la "- ll 6(d l ‘ l - \a Vd \) \x Ua - xk\^ d *'8(\a wJ \ - d“") 


(d“ l - |u lv |)! 


(\a wJ \ - d u n a )\ 


es Ik 


(A. 154) 


This is consistent with the claimed bound (4.77), and therefore finishes the proof of 
lemma 3. □ 
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